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Abstract 

We study the statistical properties of the sum St = Jq dt' a t ' , that 
is the difference of time spent positive or negative by the spin at, lo- 
cated at a given site of a Z?-dimensional Ising model evolving under 
Glauber dynamics from a random initial configuration. We investigate 
the distribution of St and the first-passage statistics (persistence) of 
this quantity. We discuss successively the three regimes of high tem- 
perature (T > T c ), criticality (T = T c ), and low temperature (T < T c ). 
We discuss in particular the question of the temperature dependence 
of the persistence exponent 6, as well as that of the spectrum of ex- 
ponents 0(x), in the low temperature phase. The probability that the 
temporal mean St/t was always larger than the equilibrium magne- 
tization is found to decay as t~ e ~2. This yields a numerical deter- 
mination of the persistence exponent 6 in the whole low temperature 
phase, in two dimensions, and above the roughening transition, in the 
low-temperature phase of the three-dimensional Ising model. 



1 Introduction 

Consider a .D-dimensional lattice of Ising spins, evolving under Glauber dy- 
namics at a fixed temperature T from a random initial configuration. This 
is obtained e.g. by quenching the system from very high temperature down 
to T. The purpose of this work is to study the influence of temperature on 
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the statistical properties of the sum 



S t = dt' a t ' , 
J o 



(1.1) 



where a t = ±1 is the spin at a given site. More precisely, we shall investigate 
two facets of this problem: 

• the scaling of St with t, and more generally the bulk properties of the 
distribution of St, 

• the statistics of rare events (such as large deviations, first passages, 
persistence), i.e. the tail properties of this distribution. 

Equivalent writings of the sum are S t = T t + — Tf — t M t , where 



are the lengths of time spent by the spin a t in the positive (negative) direc- 
tion, or occupation times of the ± states, with T t + + Tf = t, and M t is the 
temporal mean of a t , or local mean magnetization. If one views a t as the 
steps of a fictitious random walker, then S t is the position of the random 
walker at time t, and M t its mean speed. 

Simple as it may seem, the problem thus stated is actually very intricate. 
The reason is that, the values of a t at different instants of time being in 
general correlated, this study pertains to that of sums of correlated random 
variables, for which no universal result exists. The central limit theorem 
holds only for weakly correlated random variables. As we shall see, St obeys 
the central limit theorem for T > T c , while the latter is violated for T < T c . 

Hereafter, we focus our attention on the asymptotic behaviour of the 
following quantities, as t — > oo. 

(i) The probability distribution of S t , or alternatively that of M i? given by 






(1.2) 



with corresponding densities, related by 



fAt(t,x) =tf s (t,y = tx) 



d 

dx 



P(t,x). 
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Assume that the typical value of St scales, for long times, as t a . By definition 
of St, oc < 1 necessarily. Then, if a < 1, the events {S t > tx}, with x > 0, 
are rare. (Or the events {S t < tx} with x < 0.) The tail probability P(t,x), 
which measures their weight, is vanishingly small, as t — > oo. In contrast, 
if a = 1, P(t,x) is expected to converge to a limiting distribution, i.e., the 
distribution of St has fat tails. The former case corresponds to T > T c , the 
latter to T < T c , as detailed below, 
(ii) The first-passage probability 

R(t,x) =V >x, for alH' < t\ (-1 < x < 1), (1.3) 

which is the probability for the random walker not to cross the line S t > = t'x, 
up to time t. This is also referred to as the probability of persistent large 
deviations, sinceit involves the large-deviation event {S t /t > x}, and the 
persistence condition {for all t' < t}. 

Special cases of these quantities are, first, the probability of first passage 
of S t by the origin 

R(t, 0) = V{S t > > 0, for all t! < t), (1.4) 

and, secondly, the persistence probability of the process a t (assuming that 
cr t =o = +1) 

p (t) = V(no flip of o- t <, for all t! <t)=V (S t = t) , (1.5) 

which is also formally equal to R(t,x = 1). 

The principal motivation behind such an investigation comes from the 
problem of phase persistence, where, for a system undergoing phase ordering, 
the question posed is: "What is the probability for a given point of space to 
remain in the same phase as time passes?". 

For the particular case of T = 0, say for a spin model, the question 
above is answered by the knowledge of Po(t), which decays at long times 
as t~ e , thus defining the persistence exponent 9 [|TJ. However, as shown in 
[0], consideration of the more general events {St/t > x}, beyond that of 
the most extreme event {St = t} (the spin never flipped), allows both a 
stationary definition of the persistence exponent 9, as the edge singularity 
of lim^oo P(t, x), for x ±1, and the introduction of a whole spectrum 
of exponents 9(x), through the temporal decay of R(t,x). The interest of 
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introducing such concepts is strengthened by the fact that, for < T < T c , 
the same definition of the persistence exponent holds, now for x — > ±m cq 
(see below). 

In essence, the change in viewpoint when considering P(t,x) and R(t,x) 
instead of Po(t) amounts to shifting from the original question posed above to 
the more general one: "How long did the system remain in a given phase?". 
This means in particular searching the distribution of the length of time spent 
by the system inthe given phase, or occupation time of the phase, directly 
related to the sum fllTTl ). The statistics of the occupation time provides 
information on the ergodic nature of the process. Further developments on 
this theme can be found in refs. [[| to ]13| . 

This work is a sequel and a completion of 0]. We investigate the be- 
haviour of the distribution of S t and of R(t, x) in the three temperature 
regimes, T > T c , T = T c , and T < T c , and in particular revisit the question 
of the definition of the persistence exponent in the low-temperature phase. 
We will mainly consider here the two-dimensional case, and incidentally com- 
ment on the case D = 3. 



2 Temperature regimes 

In its heat-bath formulation, Glauber dynamics consists in updating the spin 
a t , located at a given site, with the probability 

PK + ^ = +l) = ^l+tanh^, (2.1) 

where h t is the local field at the given site, equal to the sum of the neighbour- 
ing spins. Under this dynamics spins thermalize in their local environment. 

The behaviour of the system is qualitatively different according to the 
value of the temperature at which the dynamics takes place, that is the 
temperature T after the quench. Three regimes are observed when letting 
T vary from high (T > T c ) to low temperatures (T < T c ). An overview 
of the changes in behaviour with temperature is given in figure [TJ, which 
depicts samples of spin histories, i.e., the position St of the Ising random 
walk as a function of time, for < t < 1000, together with the correspond- 
ing distributions of St at time t = 1000, as temperature varies. The top 
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pictures correspond to a sample of 30 spins located at various sites of a two- 
dimensional 4096 2 square lattice, the bottom pictures to all the spins of the 
lattice. (T 2D = 2/ln(l + ^2) « 2.269.) 




Figure 1: Top curves: samples of spin histories. From left to right: T = 0, 
T = 0.9 T c , T = T c , T = 1.1 T c , T = oo. The vertical axis is time (0 < 
t < 1000), the horizontal axis represents St, the position of the Ising random 
walk. Bottom curves: unnormalized histograms of St, at time t = 1000. (The 
system size is 4096 2 .) 

Snapshots at time t = 1000 of spin configurations and of configurations 
of the values of St for various temperatures are given in figures ^| and ||. 

2.1 High temperature 

As long as T > T c , thermal equilibrium is attained exponentially fast (with a 
finite relaxation time t eq ). The equilibrium magnetization m eq of the system 
is equal to 0. 

Let us first consider the simplest case T = oo. The values of the spin 
a t at successive instants of time are independent and equiprobable, as the 
steps of a binomial random walk in continuous time, hence St is a sum of 
independent, identically distributed random variables. Therefore the central 
limit theorem holds. The typical value of St reads 

(^)typ ~ *3, 
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Figure 2: Snapshots at time t = 1000 of spin configurations for various 
temperatures. From top left to bottom right: T = 0, 0.7 T c , 0.9 T c , T c , 1.1 T c , 
and T = oo. 
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Figure 3: Snapshots at time t = 1000 of configurations of the values of St 
for various temperatures. From top left to bottom right: T = 0, 0.7 T c , 
0.9 T c , T c , 1.1 T c , and T = oo. Values of S t ranging from —1000 to 1000 are 
represented by colours ranging from dark blue to dark red. Persistent spins 
(S t = ±1000), present in the top left picture (T = 0) only, are represented 
in black. 
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Figure 4: R(t = 10000, x) in the high temperature phase, for — 1 < x < 0. 
From top to bottom: T = oo, 1.5 T c , 1.3 T c , 1.1 T c , 1.05 T c , 1.01 T c . (The 
system size is 4096 2 .) 
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and the limiting distribution of St/t 1 ^ 2 is Gaussian. More precisely, the bulk 
of the distribution reads 

P(t,x = t~*z) = V (t~*S t > z) ~ j°° du exp 

because the variance of S t equals t. Hence, for (the bulk of) the density of 
M t , we have 

( tx 2 \ 

f M (t,x)~t2exv{-—\. (2.2) 

Figure [l] gives an illustration of the random walk associated to St, and of 
the Gaussian distribution of this quantity at T = oo (first pictures from the 
right). 

If t is large, and x > 0, P(t, x) measures the probability of the rare events 
such that S t ~ t, or M t ~ 0(1), i.e., of large deviations of S t with respect to 
its typical behaviour. For any sum St of independent, identically distributed 
random variables, the tail of P{t, x) (probability of a large deviation) is given 
by 

P(t,x) ~ exp(-/(x)t) (x>0). (2.3) 

t — >oo 

In the present case, T = oo, that is for a binomial random walk (with 
equiprobable steps ±1), we have I{x) = |[(l + x) ln(l + x) +(l — x) ln(l— x)}. 
(For a simple derivation see 0.) For x — > 0, I(x) ~ x 2 /2, the scaling variable 
is x 2 t = z 2 , hence the central limit theorem is recovered. 

For a sum St of independent random variables, the probability of first 
passage by the origin, V{S t > > 0, for all t' < t), decays as t - a, as is well 
known. The behaviour of R(t,x), for x ^ 0, is more subtle. According to 
ref. we have, for a binomial random walk, as t — > oo, 

!exp (-l(x)t) < x < 1 
H x = . (2.4) 

Roc(x) -1 < x < 

The asymptotic value Roo(%) is the top curve of figure f|. It is obtained by 
taking the value at time t = 10000 of R(t, x) for T = oo, depicted in figure |5] 
(left). The devil's staircase thus obtained, a discontinuous curve at all the 
rationals, is indistinguishable from its analytical prediction given in ref. |J. 
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Figure 5: R(t,x) at T = oo (left) and T = T c (right). 

For T = 00, from bottom to top: x = 1, 0.8, 0.6, 0.4, 0.2, 0.1, 0.02, 0, -0.02, 
-0.1, -0.2, -0.4, -0.6, -0.8. The dashed line parallel to R(t,Q) has slope 
-0(0) = -§. 

For T = T c , from bottom to top: x =1, 0.8, 0.6, 0.4, 0.2, 0.1, 0, -0.1, -0.2, 
—0.4, —0.6, —0.8. The dashed line parallel to R(t, 0) has slope — 8(0) = 
-0.273. (The system size is 4096 2 .) 
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The persistence probability pa(t) behaves, as t — > oo, as po(t) ~ e -( ln2 ) t ) 
which matches both Q2.3|) and ( |2.4|) for x = 1. Note that (|2.4j ) defines the 
first-passage exponent 6(x = 0) = ~. 

Let us now investigate the deformations induced on the quantities con- 
sidered above, when T decreases. 

If T c < T < oo, the two-time autocorrelation (a s a t ) (with s < t) is short 
ranged. At equilibrium, i.e., for 1 C s, (o~ s cr t ) ~ exp (— (i — s)/t eq ). Hence 
the statistics of St is that of a random walk with short-ranged correlations, 
and the central limit theorem still holds. The distribution of St is Gaussian 
(see second pictures from the right in figure 1, for T = 1.1 T c ), with a width 
proportional to (it e q) 5 - While the Gaussian shape of P(t,x) in the bulk is 
universal, as long as the correlations (<7 s cr t ) are short-ranged, the tails are 
not. Hence I(x) is expected to be deformed as T decreases from infinity to 
T 

^ c- 

The probability of first passage by the origin, V(S t > > 0, for all t' < t), 
still decays as t~2. The persistence probability R(t,x) behaves qualitatively 
as in equation ( |2.4|) . As T decreases down to T c , the discontinuous curve 
-Roo(x) is deformed, as seen in figure £|. Gaps vanish when T — > T c . Note the 
finite time effect when approaching T c , manifested by the fact that Roo(x) 
does not vanish as x — ► 0. 

2.2 Critical coarsening 

The system is now quenched from a disordered initial state to its critical 
point. 

After the quench, spatial correlations develop in the system, just as in 
the critical state, but only over a length scale which grows like t 1 ^ Zc , where 
z c is the dynamic critical exponent. The equal-time correlation function has 
the scaling form 

C x (t) = (a t (O) ( T i (x)) = |x|- 2 ^0 

where x is a lattice point, and /3 and v are the usual static critical exponents. 
The scaling function <p(y) goes to a constant for y — > 0, while it falls off 
exponentially to zero for x — > oo, i.e., on scales smaller than t 1//Zc the system 
looks critical, while on larger scales it is disordered. In two dimensions z c w 
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2.17, j3 — |, and v = 1; in three dimensions 2 C ~ 2.04, /3 ps 0.327, and 
i/ « 0.63 @. 



2 




Figure 6: Critical scaling function <£> obtained by rescaling (t, x) for three 
different times: t =100, 1000, 10000, see equation (|2.7|) . (The system size is 
4096 2 .) 

For 1 C s ~ i the two-time autocorrelation function scales as [|l|, [16[ 

C(t, s) = (a s a t ) = s- 2 ^g c (~) . (2.5) 



A numerical determination of the scaling function g c can be found in ref. fllC 
This implies, setting a = 2/3/ z/2 c , 



Mr 



2 /•' 



ft rts 

/ dt 2 / cZtj (cr tl o- t2 ) 
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= 2t~ a jf dv £ du u~ a g c (^j 
= r a [j^£dxx- a g c {x)y 

hence 

(^)typ ~ t 1 -^. (2.6) 

Note that this scaling behaviour was obtained from the non- stationary form 
( |2.5| ). The behaviour of S t at criticality is illustrated by the middle pictures 
of figure 1. 

From equation ( |2.6| ) one can infer that the bulk of the probability density 
of M t scales as 

f M {t 1 x)~t?/"*< v ,(tV 1 "°x), (2.7) 

which is the critical counterpart of ( |2.2| ). The scaling function (p is depicted 
in figure |6|, for D = 2, with (3/vz c « 0.0576. Due to the smallness of this 
exponent, ju is very slowly peaking. The existence of the scaling form ( |2.7j ) 
implies gap scaling for higher moments, i.e., ^M| fc \ ~ (M t 2 ) fc . 

The probability R(t, x) behaves qualitatively as in the high temperature 
phase (see figure |, right). For < x < 1, it decreases to zero faster than a 
power-law (in particular po(t) = R(t,x = 1)). For — 1 < x < 0, it decreases, 
extremely slowly, to a constant Roo(x) obtained by extrapolating the results 
of figure|, asT -> T c . Finally for x = 0, R(t,x) ~ t~ e ^. We find 9(0) « 0.27 
in 2D and 6(0) ~ 0.41 in 3D. Note that ^(0) is a nonequilibrium critical 
exponent. 

Remark. The exponent #(0) can also be measured at criticality for the 
two-dimensional pi —p2 model |17|, |j"8f . This model is defined as follows. The 
right side of (£J), p(h) = |(l+tanh/i/r), is a function of h (h = ±4, ±2,0), 
such that p(—h) = 1 — p(h). In particular p(0) = |. The dynamics of the 
Ising model therefore depends on one parameter only, which is p\ = p(2), or 
alternatively p 2 = p(4), the two being related. Considering instead these two 
quantities as independent (\ < pi,p 2 < 1) defines the p\ — pi model. For 
the voter model (p\ = |, P2 = 1) one finds #(0) ~ 0.37. As one moves along 
the critical line, from thevoter point to the Ising critical point (p± = 0.854, 
P2 = 0.971), numerical measurements show algebraic decay of R(t,0), with 
a seemingly constant slope on a log-log plot (i.e. no sign of a crossover 
to 9(0) ~ 0.27), defining a decreasingexponent #(0). Then, from the Ising 
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critical point to the end of the curve (pi = 1, P2 ~ 0.85), this exponent is 
constant, 9(0) « 0.27. 

2.3 Low temperature coarsening 

In the low-temperature phase (0 < T < T c ), the lower the temperature, 
the higher the tendency of a spin to align with the majority. Hence the 
system coarsens, i.e., domains of opposite signs grow, because the system 
tries to reach locally one of the two equilibrium phases, corresponding to an 
equilibrium magnetization ±m eq , where (in 2D) 

m eq (T) = (l-(sinh2/T)- 

In the scaling regime the system is statistically self-similar, with only one 
single characteristic length scale, which is the size of a typical domain. 

For a typical spin, deep into a domain, two scales of time are observed: a 
fast one, due to thermal flips, and a slow one, due to the passage of domain 
walls. The net result of the thermal flips can be seen on figure 1. For 
example, for T = 0.9 T c (second pictures from the left), the distribution of St 
concentrates, i.e., no longer covers the whole interval [—£,£]. A more precise 
statement is given presently. 

In the regime 1 C s ~ t, the autocorrelation function scales as 

C(t, s) = r< (T) g . (2.8) 



According to this hypothesis [|19j , all the temperature dependence is factored 
out in the prefactor m^. Thus, by the same computation as above, 



(M 2 ) = m e 2 q (T) j Q dxg(x). (2.9) 

Hence 

(St)typ ~ *, 

which indicates the existence of a limiting distribution for M t , as t — > oo: 

Im(x) = lim f M {t,x). 

t — >oo 
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This distribution is observed to be a U-shaped curve defined on [— m eq ,m eqJ 



Equation ( |2.9| ) shows that the variance of M t /m eq no longer depends on 
temperature. Assuming that the same factorization of m eq holds for higher 
moments implies that the whole limiting distribution of the rescaled variable 
M t /m eq is independent of temperature, and therefore identical to the T = 
distribution. The latter being singular at x — > ±1, with singularity exponent 
equal to 9 — 1 0, [|, the existence of a unique master curve for the distribution 
of M/m eq thus provides a simple definition of the persistence exponent 9 at 
finite temperature. 

As seen in figure ^ the zero-temperature limiting distribution ju is 
extremely close to a beta distribution ~ (1 — x 2 ) 0-1 , hence, in the low- 
temperature phase, we have 



x 2 



In two dimensions, the determination of 9 from fu at T = yields 9 ~ 0.22 
(see figure |7|), confirming previous estimates obtained either by numerical 
simulations [201 , or in an experiment on a liquid crystal system (See 
also HO 

Parallel observations come from the measurement of _R(t, x). For — m eq < 
x < m oq , we have 



i?(t, ar) = P(S t ,/t' > x for all t' < t) ~ t 



-«(*) 



which defines the spectrum of exponents #(x) 0, [J. Outside this interval, 
i?(t, x) either decays faster than a power-law, or goes to a constant (see 
figure [8p. In two dimensions we observed with reasonable accuracy that 
9(x/m eq ) does not depend on temperature. This generalizes the observa- 
tion, made in 0, for the two-dimensional Ising model, that the first-passage 
exponent #(0) (~ 0.19 in 2D) is independent of temperature for any T < T c . 
The case x = m eq deserves special mention. We find 

R(t,m cq ) = V(S t ,/t' > m cq for all t' <t) ~ t~°~^, 
which is equivalent to saying that 

6{m eq )=e + - (0<T<T c ). 
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Figure 7: Probability density f M (t,x) at T = for t = 100, 1000, 10000, 
20000. Inset: log-log plot of the same curves against 1 — x 2 . The dashed line 
has slope —0.78. (The system size is 4096 2 .) 
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Figure 8: R(t,x) at T = 0.895 T c (corresponding to m eq = 0.9) (left) and 
T = (right). 

For T = 0.895 T c , from bottom to top: x = 1, 0.96, 0.94, 0.92, 0.9, 0.88, 0.82, 
0.5, 0, -0.6, -0.9. The dashed line has slope -6(m cq ) = -0.22 - \. 
For T = 0, from bottom to top: x = 1, 0.8, 0.6, 0.4, 0.2, 0, -0.2, -0.4, -0.6, 
-0.8, -0.9. (The system size is 4096 2 .) 
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Figure 9: R(t,x) at T = 0.75 T c (corresponding to m eq ~ 0.9) for the three- 
dimensional Ising model. From bottom to top: x — 1, 0.9, 0.8, 0.6, 0.4, 
0.2, 0, -0.2, -0.4, -0.6, -0.8, -0.9. The dashed line has slope -0(m eq ) = 
—0.26 — |. (The system size is 512 3 .) 
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The presence of the exponent | can be simply understood. The decay of 

R(t, m eq ), i.e., of the probability for the Ising random walker not to cross the 

line S t = m eq t has two causes. On the longer scale it is due to the crossing of 

domain walls, on the shorter time scale it is due to thermal fluctuations of M t 

around m eq . These fluctuations are the same as at equilibrium. For a system 

at equilibrium the random walk is biased, since in average (M t ) = m eq . The 

probability for this random walk not to cross the line S t = m eq t decays as 
i 

t-2. 

Illustrations of this phenomenon are given in figure [8] for the two-dimen- 
sional case, and in figure |9] for a three-dimensional cubic lattice. In these 
figures R(t, m eq ) appears as a separatrix, with slope — 9 — |, where 9 ~ 0.22 in 
2D, and 9 ~ 0.26 in 3D. This last value is in agreement with the measurement 



of refs. [[23|, p4j| . Note that the accuracy is better in 3D than in 2D, though 



the size used in the former case is much smaller. The reason is that, in 
order to observe the factor i~2, m eq has to be significantly different from 1, 
and therefore T close enough to T c . On the other hand, if T is too close 
to T c , crossover effects become important and preclude the observation of 
algebraicity of R(t,x). Due to the temperature dependence of m eq , a better 
compromise is found in 3D. For instance m eq = 0.9 for T = 0.895 T c in 2D, 
while this value is already obtained for T = 0.75 T c in 3D. (T C 3D m 4.511 

[00 

Note that the temperature chosen, T = 0.75 T c , is above the roughening 



transition of the 3D Ising model, Tr ~ 0.544 T c p6| . Below Tr the exponents 
9(x), and in particular 9, take smaller values. 

The persistence probability po(i) is algebraically decaying as t~ e at T = 
only. 



3 Discussion 

Let us give a brief summary. The three temperature regimes, T > T c , T = T c 
and T < T c , correspond, for the process a t , to correlations of increasing 
strength, and, correspondingly, to different behaviours of thesum St. 

• At high temperature, for 1 < s ~ the two-time autocorrelation 
function C(t, s) = (<J s <7 t ) (s < t) is short ranged and stationary: 



C(t,s) ~ e 
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hence St obeys the central limit theorem, and 



(^)typ~^. 



• At criticality, for 1 C s ~ t, 



C(t,s) = s 



2j3/vz c 




) 



implying 




~ t 



1-P/VZC 



• In the low-temperature phase, for 1 <C s ~ t, 



C(t,s)=ml q (T)g^. 



implying 




~ t. 



In the first two cases (T > T c ), the probability density of S t /t = M t 
is peaking, while that of the scaling variable St/(St)typ is converging to a 
Gaussian, for T > T c , and to the function <p (cf eq. ( |2.7| )), for T = T c . 

Below T c , the probability density /m(A x) concentrates on [— m eq , m eq ], as 
t — > oo. In the variable x/m eq , it can be rescaled onto a universal curve, which 
is the T = limiting distribution Jm, singular at x — > ±1, with singularity 
exponent equal to — 1 (see figure [7]). Another alternative definition of 
the persistence exponent at finite temperature is provided by the decay of 
R(t,m eq ) ~ t~ e ~\. 

In summary, the central limit theorem for St is violated for any T < T c . 
On the other hand ergodicity is broken for T < T c , since M t , the temporal 
mean of the spin at, remains distributed, as t — > oo, instead of converging to 
the average (a t ). 

We conclude by a last comment. The central limit theorem can be violated 
in essentially two ways. Either by adding identically distributed random 
variables, all of the same orderof magnitude, i.e. with a narrow common 
distribution, but otherwise strongly correlated. Or by adding independent 
identically distributed random variables, with a broad common distribution. 
The violation of the central limit theorem for St at T = can be viewed in 
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either of these two ways. On one hand, by its very definition, St is the sum of 
the correlated random variables at, as discussed in this paper. On the other 
hand, St can also be viewed as the alternating sum of the intervals of time 
Ti, r 2 ..., between two sign changes of a t , which are broadly distributed. 
Though these intervals of time are not independent, the violation of the 
central limit theorem is nevertheless essentially due to the algebraic tail ~ 
r -( 1 + 6 f their distribution. The case < T < T c is more difficult to interpret 
this way because of the occurrence of thermal flips. A natural model to 
consider however, is one where the intervals of time t, are independent with 
distribution p{r) ~ r~ l ~ e ||. Though this is a simplification of reality, it 
leads nevertheless to the same classification of behaviours for C(t,s) and 
(<St)typ as summarized above, with respectively 6 > 2 corresponding to T > 
T c , 1 < 9 < 2 to T = T c and 6 < 1 to T < T c [|Kj. 
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